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ABSTRACT 

We present new X-ray timing and spectral results on the 8.0-second X-ray pulsar CXOU J010043.1- 
721134 from a series of observations using the Chandra X-ray Observatory. We find a spin period 
in 2004 January of 8.020392 ± 0.000009 seconds. Comparison of this to 2001 Chandra observations 
implies a period derivative P — (1.88 ± 0.08) x 10~ n s s^ 1 , leading to an inferred dipole surface 
magnetic field of 3.9 x 10 14 G. The spectrum is well fit to an absorbed blackbody of temperature 
kT = 0.38 ± 0.02 keV with a power law tail of photon index T — 2.0 ± 0.6. We find that the source 
has an unabsorbed X-ray flux (0.5-10 keV) of 4+ 2 x 10~ 13 erg cm~ 2 s _1 and a corresponding X-ray 
luminosity of ~ 2 x 10 35 erg s _1 for a distance of 60 kpc. These properties support classification of 
CXOU J010043. 1-721134 as the seventh confirmed anomalous X-ray pulsar, the eleventh confirmed 
magnetar, and the first magnetar to be identified in the Small Magellanic Cloud. 
Subject headings: pulsars: general — pulsars: individual(CXOU J010043. 1-721134) — stars: neutron 



1. INTRODUCTION 

Before the mid-1990s, it was thought that young neu- 
tron stars had short (P <C Is) spin periods, surface mag- 
netic fields B~ 10 12 G, and exhibited radio pulsations. 
However, recent theoretical and observational work has 
revised this picture. We now know that a small popula- 
tion of sources known as soft gamma repeaters (SGRs) 
are young neutron stars with extreme (B~ 10 14 ~ 15 G) 
surface magnetic fields, or "magnetars". Another class 
of objects, anomalous X-ray pulsars (AXPs), are thought 
to be closely related to SGRs. Anomalous X-ray pul- 
sars are so-named due to their high X-ray luminosi- 
ties and unusually fast spindown rates, with no evi- 
dence of variation due to binary motion, which distin- 
guish them from both isolated radio pulsars and accret- 
ing X-ray binaries. Most significantly, AXPs exhibit a 
small range of spin periods (5-12 seconds) with steady 
spindown that, if powered by magnetic dipole radia- 
tion, imply magnetic fields 2-3 ord ers of magnitude larger 
than those of typical radio pulsars ijThompson fc Duncan! 
1996). The magnetar model for AXPs is supported by 
the SGR-like bursting activit y that has been observed 
from some of these so urces l|Gavriil. Kasni. fc Woods! 
I2002t iKasni et alJ 120031) . Other characteristics include 
spectra that are well fit by a blackbody model with a 
temperature kT ~ 0.5 keV and a power law tail with 
photon index T=2-4, and lack of a measurable binary 
companion that could otherwise p roduce similar X-ray 
lumin osities through accretion fsee lWoods fc Thompson! 
20051 for a summary of the properties of known mag- 
netars). With unabsorbed luminosities in the range of 
10 34 to 10 36 erg s _1 , inferred magnetic fields of 10 14 - 
10 15 Gauss, and characteristic ages of ~ 8 — 200 kyr, 
it has been firmly established that AXPs belong to the 
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magnetar family ijOavriil. Kasoi. fc Woodsll2004D . 

The X-ray source CXOU J010043.1-721134 in the 
Small Magellanic Cloud (SMC) was first detected by 
Einstein in 1979, and was later observed several times 
by ROSAT and ASCA, but its periodicity was not no- 
ticed until an archival search for X-ray pulsations was 
carri ed out on a 2001 C handra observation of a nearby 
field l|Lamb et al] 120021) . The period of this source was 
determined to be 8.0 seconds and its X-ray luminosity 
was measured to be on the order of 10 35 erg s" 1 , assum- 
mg a distance to the SMC of 60 kpc jLamb et al.1 '2002. 
120031 iMaiid. Lamb, fc Macombll2004|) . Such a high lu- 
minosity and long spin period suggested tha t this object 
may belong to the growing class of A XPs ijLamb et alJ 
20021 IMaiid. Lamb, fc Macomrll27)04j) . 

In this Letter we present new Chandra results support- 
ing CXOU J010043. 1-721134 as the newest member of 
the AXP class. We present timing analysis from which 
we measure a period derivative, spectral results, and an 
analysis of archival Chandra and XMM data that as- 
sesses the flux characteristics of the source over time. 

2. OBSERVATIONAL SETUP 

Five new Chandra observations of CXOU J010043.1- 
721134 were completed between 2004 January 27 and 
2004 March 25. The source was clearly detected at each 
epoch, located at (J2000) right ascension 01 ?i 00 m 43!03, 
declination — 72°ll'33'/6, with a radial uncertainty of 
0'.'5. Coordinates from the Chandra pointings have not 
yet been refined using an optical reference. Each observa- 
tion had approximately 16 ks of on-source time, and the 
observations were spaced in a geometric series designed 
for a phase-coherent timing analysis. All five observa- 
tions were made with the ACIS-S3 CCD in Very Faint 
mode. A 1/8 subarray was used to reduce the time res- 
olution to 0.4 s. 

Additionally, we analyzed serendipitous archival ob- 
servations from both Chandra and XMM. The archival 
2001 May 15 Chandra observation was made in faint 
mode, and the source fell on the ACIS-I1 CCD approxi- 
mately 10' from the aimpoint. We reprocessed the level 1 
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events file and applied the CTI, standard gain, and time- 
dependent gain corrections. XMM archival observations 
were taken on 2000 October 10 and on 2001 November 
20 using the EPIC MOS detectors, and the 2001 obser- 
vation also used the EPIC PN detector. Both datasets 
lost some observation time due to high background. The 
source location in each XMM observations is also off- 
axis, and the source lies on a chip gap on the MOS2 
instrument. Therefore, we discarded all MOS2 data in 
this analysis. Table ^ contains detailed information on 
each observation. We processed Chandra observations 
using CIAO 3.1 and CALDB 2.28, and XMM observa- 
tions using SAS 6.0.0. 

3. TIMING ANALYSIS 

The five 2004 Chandra observations of 
CXOU J010043. 1-721134 were separated by suc- 
cessive intervals of 0.40, 1.72, 8.9, and 46.4 days. 
For each epoch, we extracted the events from a 2" 
radius region centered on the coordinates of the source 
and corrected their arrival times to the solar system 
barycenter. We measured pulse profiles for each ob- 
servation by folding the X-ray events at a range of 
periods around the P~8. 02s spin-per iod of the pulsar as 
reported bv lLamb et all l|2002. 12003). The profiles with 
the highest significance (as determined by maximizing 
X 2 with respect to a model with no pulsations) were 
cross-correlated and co-added in phase to create an 
integrated template pulse profile from all the obser- 
vations. We then determined pulse arrival times by 
cross-correlating the profiles from each observation with 
the high signal-to-noise template profile. From the three 
Chandra observations in 2004 January, we unambigu- 
ously measure a spin period P = 8.020392(9) s, where 
the number in parentheses indicates the uncertainty in 
the last digit. Unfortunately, the fractional phase errors 
from these measurements were ~0.06, significantly larger 
than the 0.02—0.03 fractional phase errors that we were 
expecting. These larger arrival time uncertainties mean 
that the spin period determined above from the three 
2004 January observations was not of sufficient accuracy 
to unambiguously account for each rotation of the pulsar 
between the 2004 January and February observations, 
let alone during the much longer gap between the 2004 
February and March observations. 

From the archival Chandra observation taken on 2001 
May 15, we measure P — 8.0188(1) s, which is signifi- 
cantly shorter than P = 8.020392(9) s in 2004 January as 
determined above. If this period difference is due to a rel- 
atively constant spin-down, it implies P=1.9(l)xl0~ n 
s s _1 . If we then force a phasing of the five 2004 
observations to most closely match the average spin- 
down between 2001 May and 2004 January, we determine 
P=1.88(8) x 10~ n s s _1 . A single phase wrap between 
the 2004 February and 2004 March observations implies 
P=1.26(8)xl(T u s s" 1 or P=2.50(8) x 10" 11 s s -1 , both 
of which are clearly inconsistent with the average spin- 
down. It is therefore highly likely that the true spin-down 
value of this source is P=1.88(8) x 10" 11 s s" 1 . 

4. SPECTROSCOPY 

For the spectral analysis on the 2004 Chandra obser- 
vations, we used an extraction radius of 2" for the source 
region and 29'.'03 for the background region. For the 



2001 off-axis Chandra observation, the source region was 
an ellipse with major and minor axes 10'.'82 and 9'.'84 re- 
spectively, and background extraction circle with radius 
92'.'50. In all Chandra analyses, the background regions 
were not centered on the source but did cover the source 
and so the source region was excluded. For the XMM 
observations, the source extraction circles for the MOS 
and PN detectors had radii of 12'.' 40 and 19'.'20, respec- 
tively, while the background regions for both were 65'.'60 
and did not overlap the source region. 

We successfully fit each spectrum in Sherpa 5 to a 
model containing an absorbed blackbody plus a power 
law, accounting for absorption from the Milky Way and 
the SMC. We assumed a Galactic absor bing column 
N H =4.3 x 10 20 cm" 2 (jNaze et alJl2003|) . Standard 
SMC elemental abundances were used the hydrogen col- 
umn due to the SMC co ntribution was allowed to vary 6 
ijRussell fc Dopitai rr992'). The temperature and photon 
index, and their normalizations, were allowed to vary in 
their respective models. 

Table [21 compares various models that we fit to the 
data. A blackbody plus power law model provides a 
much better fit tha n the single model fits that were pre- 
viously attempted lEamb et al l 12001 iNaze et alJ 120031: 
iMaiid. Lamb, k. Macombll2004|) . Fitting to a simple ab- 
sorbed power law model yields a poor result, as shown 
by the reduced x 2 = 2.12. Fitting to a simple absorbed 
blackbody model yields a reduced x 2 = 1-10- How- 
ever, this fit is infeasible because it requires that there 
be no absorption column contribution from the SMC 
(Nh smc)- I n fact, if Nhsmc i s frozen at a reasonable 
value we are still not able to fit a sensible model to the 
data. 

Statistics were calculated to 90% confidence (<r = 1.6). 
In each dataset, we excluded data below 0.5 keV to re- 
flect telescope sensitivity and above 4.0 keV because the 
source spectrum is not detected beyond this. We sub- 
tracted the background, although it was minimal in the 
Chandra data, and grouped the spectra into bins of at 
least 25 counts for analysis. 

4.1. Chandra 2004 Observations 

We fit the five 2004 Chandra datasets simultaneously 
and find that N HSMC = (3 ± 4) X 10 21 cm" 2 . The 
large uncertainty in this absorbing column is due to the 
coupling between the absorbing column and the power 
law index T. We find a blackbody temperature kT = 
0.38±0.02 keV, power law photon index T = 2.0±0.6, 
and reduced x 2 — 0.97. The unabsorbed 0.5-10 keV lu- 
minosity of the source in these Chandra observations is 
~ 2 x 10 35 erg s _1 at a distance of 60 kpc, and the black- 
body component contributes ^50% of the luminosity. 

4.2. 2000 and 2001 Archival Observation 

We also analyzed the 2001 archival Chandra obser- 
vation. The model applied to the previously described 

5 http://cxc.harvard.edu/sherpa/ 

6 Abu ndances for elemen t s Co a nd Ca were not included in the 
model of Russell & DoDita (1992), and so both were assigned a 
value of 0.3171 with respect to solar in the spectral model based 
on the average abundance of the other elements. Variations in the 
Co and Ca abundances do not result in significant changes to the 
fits. 
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Chandra data also yields an excellent fit to this observa- 
tion, with Nhsmc = (8±6)xl0 21 cm" 2 , kT = 0.36±0.03 
keV, T = 1.9±0.4, and reduced x 2 = 1-15. The unab- 
sorbed 0.5-10 keV luminosity of the source is ~ 2 x 10 35 
erg s , and the blackbody component contributes ^60% 
of the luminosity. It should be noted that these data were 
fit after applying a user-contributed time-dependent gain 
correction. This software accounts for an instrumental 
effect that has not been fully addressed by the default 
Sherpa package for older observations. 

The best fitting model to the 2000 and 2001 XMM 
data is again the absorbed blackbody with a power law 
component. Though we also fit the two XMM datasets 
individually, we found that they were well modeled when 
fit jointly. This model implies N H smc = (5 ± 7) x 10 21 
cm" 2 , kT = 0.31±0.04 keV, T = 2.0±0.5, and reduced 
X 2 = 1.08. The unabsorbed 0.5-10 keV luminosity of the 
source as seen by XMM is ~ 2 x 10 35 erg s _1 , and the 
blackbody component contributes ~40% of the luminos- 
ity. 

5. LONG TERM FLUX BEHAVIOR 

We have also examined the flux variability of 
CXOU J010043. 1-721134 over the last three decades us- 
ing the dat a in j^l plus count rates from archival data 
reported bv lLamb et al.l l|2002j) . We inferred the fluxes 
and their errors for the older observations using the fol- 
lowing method. For both the Chandra and XMM obser- 
vations presented in this paper, we used Sherpa's eflux 
command to calculate the flux for the best-fit parameters 
of Nhsmc^ r, and kT. We also determined minimum and 
maximum possible fluxes, by applying the errors in the 
spectral fits, which were calculated at the 90% confidence 
level, in order to achieve the largest possible flux range. 
Specifically, the minimum flux occurs when Nh smc an d 
F are at a maximum and k T is at a minimum, and the 
maximum flux occurs in the reverse scenario. We ap- 
plied a correction to account for the correlations in the 
uncertainties of the three parameters and defined the flux 
error to be the difference between the flux defined by the 
best-fit parameters and these extreme flux values. 

Fluxes for the observations fr om previous X-ray satel- 
lites listed in lLamb et al.1 l)2002t) could not be calculated 
in this manner without reprocessing each dataset, so we 
developed a method by which to extrapolate their fluxes 
and errors. We converted the errors from our Chandra 
observations into fractional uncertainties and then com- 
pared them to the fractional uncertainty in their mea- 
sured count rate to obtain a count rate-flux uncertainty 
relation. Using the best-fit model for the 2004 Chandra 
observations (see Tabled, we converted from count rate 
to flux using HE AS ARC's WebPIMMS 7 . We then used 
our count rate-flux relation to extrapolate the uncertain- 
ties in these older observations. Results are plotted in 
Figure ^ Because small changes in the fit parameters 
lead to large variation in flux for a given count rate, it is 
difficult to compare fluxes across multiple epochs. How- 
ever, at face value, the flux behavior does not seem to 
demonstrate variation spanning more than one order of 
magnitude over a timescale of ~25 years. 

6. DISCUSSION 

7 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 



The properties of CXOU J010043. 1-721 134 are com- 
pletely consistent with those expected of an AXP. In 2004 
January, the source had a period of 8.020392(9) seconds 
and a P = 1.88(8) x 10~ n s s -1 . If this spin-down is due 
to "standard" dipole radiation, the implied surface mag- 
netic field strength is 3.2x 10 19 (PP) 1 / 2 G = 3.9xl0 14 G, 
consistent with magnetic fields of other AXPs. Using the 
estimated period derivative, it is possible to calculate the 
characteristic age and spindown luminosity of the source. 
A period derivative P — 1.88(8) x 10~ n s s _1 implies 
a characteristic age of P/2P = 6800 years (with a 4% 
uncert ainty), consistent with characteristic ages of mag- 
netars l|Gaensler et al.ll200l|) . CXOU J010043. 1-721 134 
has a spin-down luminosity of 4ir 2 IP '/P 3, = 1 x 10 33 erg 
s _1 , where I = 10 45 g cm 2 is the assumed moment of 
inertia of the neutron star. 

The source spectrum is well described by a photon in- 
dex of r = 2.0±0.6 and a blackbody temperature of kT 
= 0.38±0.02 keV. This implies a 0.5-10 keV unabsorbed 
luminosity of ^2 x 10 35 erg s _1 , where the blackbody con- 
tributes ^50% of the luminosity (at a distance of 60 kpc). 
This fit is consistent with those of other AXPs, which 
tend to have kT ~ 0.5 keV and T ~2-4. Although the 
blackbody contribution to the luminosity of an AXP is 
highly energy dependent, the value of ~50% in the case of 
CXOU J010043. 1-72 1134 is consistent with other AXPs 
in the 1-4 keV range (jOzel. Psaltis. fc Kaspll200H) . 

It has been suggested that the source could be a Be 
X-ray binary with a ~25 day period b ecause there is a n 
optical Be star within 2" of the source (N aze et all2 003 h 
X-ray binaries are also described by a blackbody plus 
power law model, but generally their blackbody temper - 
ature is higher (i.e. 1-2 keV) l|Haberl fc Pietschll2004[) . 
However, our timing analysis rules out a 10-30 day binary 
due to the absence of acceleration by a massive compan- 
ion in the first three Chandra obsevations since the phase 
connects unambiguously over ^3 days. The phase wrap 
after the third observation does leave open the possibil- 
ity of longer period binary motion, however, the long 
term period change seems to match the best phase con- 
nected solution, which likely rules out any binary com- 
panion and argues strongly for a relatively steady spin- 
down. We further note that iDurant fc van Kerkwiild 
(2005) have identified a likely optical counterpart to 
CXOU J010043.1-721134 distinct from the adjacent B 
star. This object provides a convincing alternative to 
the Be X-ray binary interpretation. 

CXOU J010043.1-721134 does not show strong evi- 
dence of variability over the timescale of years. Al- 
though not all AXPs are variable, some AXPs have been 
found to vary across one to two orders of magnitude 
l|Gavriil. Kasoi. fc Woodll2004t librahim et al.ll2004fl . X- 
ray binaries also exhibit variations, by factors of up to 20 
over a similar timescale, and can cha nge by a factor of 
more than 100 during giant outbursts l|Haberl fc Pietschl 
2004). Even considering the large uncertainties, we do 
not see variability on this scale in CXOU J010043.1- 
721134 over the past 25 years. While it does not con- 
clusively characterize the object, apparent lack of vari- 
ability across orders of magnitude does not challenge the 
categorization of this source as an AXP. 

In summary, the spectral, timing, and photometric 
properties of CXOU J010043. 1-721134 are consistent 
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TABLE 1 

Spectral Fits to Chandra and XMM Datasets 
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2.0(3) 


0.30(5) 
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the amount of usable time 


recovered from the original 
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caused by CCD readout, 
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default energy range for each instrument (Chandra: 0.2-10keV, XMM: 0.4-10keV); JV H g M £ = the absorbing column due to the Small Magellanic 
Cloud only, assuming a Galactic absorbing column of 4.3 X 10 20 cm -2 ; V = the power law photon index; kT BB = blackbody characteristic 
temperature; x 2 value is for the fit of the absorbed black body plus power law model; The unabsorbed 0.5-10keV flux for the five 2004 Chandra 
observations is 4+ 2 X 10~ 13 erg cm -2 s _1 , and in all other observations is 5 + 2 X 10~ 13 erg cm -2 s — 

a Values for MOS1 and PN respectively. 



with the properties of other AXPs l)Mereghetti et all 
2002), expanding this class of magnetars to seven con- 
firmed members. We note that the source is the only 
known magnetar in the SMC. If t he magnetar birth 
rate i s ~10% of that of radi o pulsars ijKouveliotou et aD 
179941 ICaensler et aTJ 12005). then for an SMC super- 
nova rate of (5 ± 3) x 10~ 4 per year l)Crawford et alJ 
I2001D and a magnetar lifetime on the order of 10 4 yr 
l|Colpi. Geppert. k Paee1l200Cit iGaensler et al]l200"H) . we 
can expect 0.5 ± 0.3 magnetars to exist in the SMC 
at a given time. This would be consistent with the 
fact that only this magnetar is known in the SMC thus 
far, despite several archival searches for X-ray sources 
llMaiid. Lamb k Macomr3l200l iHaberl & PietscipOOl 
lLavcock et al.ll2005j) . Interestingly. l|Ibrahim et al.ll2004(i 
report the discovery of a transient magnetar, detected 
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only because of an outburst that resulted in a substantial 
flux enhancement that faded on a time scale of months. 
Similarly, the highly variable 7-s pulsar AX J1845-0258 
may be another such source (Vasisht et al. 2000). These 
sources suggest that the magnetar birthrate could be 
substantially higher than has been estimated, possibly 
even comparable to the radio pulsar birthrate. Contin- 
ued monitoring of the SMC will provide a useful test 
of these possibilities. In any case, CXOU J010043.1- 
721134, with a well-known distance and low extinction, 
will be a valuable probe into this exotic class of objects. 
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TABLE 2 

Joint Spectral Fits to 2004 Chandra Data 
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235/214 = 


1.1 


PL 


8(l)xl0 21 


2.20(5)/ ■■■ 


~6 


455/214 = 


2.1 


BB+PL 


3(4) xlO 21 


2.0(6)/0.38(2) 


~4 


205/212 = 


0.97 



Note. — Parentheses following each value represent the error in the last digit; Njj — 
the absorbing column due to the Small Magellanic Cloud only, assuming a Galactic absorbing 
column of 4.3 X 10 20 cm -2 ; T — the power law photon index; kTsB — blackbody characteristic 
temperature; F — Unabsorbcd flux for the model in the 0.5-10 keV range. 



a The data could not be fit within Sherpa's limits for Nhsmc- 
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Fig. 1. — Estimated unabsorbed fluxes for CXOU J010043. 1-721134 over the last three decades of observations, using the BB+PL model 
from Tabic [5] in the energy range 0.2-10 keV. Each instrument is denoted with a separate plot symbol as follows: X - Einstein IPC, o - 
ROSAT PSPC, A - ROSAT HRI, □ - ASCA CIS, * - XMM MOS/PN, • - Chandra ACIS. Flux estimates from observations presented 
in this Letter were determ i ned b y using Sherpa's eflux command. Those not presented in this Letter are derived from the observed count 
rates listed in lLamb etaTI £2002) by using HEASARC WebPIMMS. 



